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There is substantial interest in producing sam- 
ples of ultracold molecules for possible applica- 
tions in quantum computation [l], 0), quantum 
simulation of condensed matter systems 0, Q], 
precision measurements 0], controlled chem- 
istry [3, @], and high precision spectroscopy [9]. 
A crucial step to obtaining large samples of ultra- 
cold, trapped molecules is developing a means to 
bridge the gap between typical molecular source 
velocities (~ 150 — 600^) and velocities for which 
trap loading or confinement is possible (< 5—20^). 
Here we show deceleration of a beam of neutral 
strontium monofluoride (SrF) molecules using ra- 
diative force. Under certain conditions, the decel- 
eration results in a substantial flux of molecules 
with velocities < 50-^^. The observed slowing, 
from ^140^, corresponds to scattering >10"' pho- 
tons. We also observe longitudinal velocity com- 
pression under different conditions. Co mb ined 
with molecular laser cooling techniques 
this lays the groundwork to create slow and cold 
molecular beams suitable for trap loading. 

Tremendous advances have been made in the decelera- 
tion of molecular beams in the past decade. Stark decel- 
eration [13 [^eeman deceleration [3l , and collisional de- 
celeration 15[ have all been demonstrated to slow molec- 
ular beams. However, only for fairly light species 20 
amu) with substantial vapor pressure at room temper- 
ature have these methods been demonstrated to allow 
slowing to velocities as low as a few m/s, as needed to 
make trapping possible (T^j. Optical deceleration has 
been demonstrated to slow molecular beams to rest [111 , 
but the high laser intensities required limit application 
to small volumes. While these methods are useful, all of 
them conserve phase-space density and hence slow with- 
out cooling. However, recently it has been demonstrated 
that radiative forces can be used to cool molecules [l^ . 
Assuming a given species is amenable to laser cooling, the 
same radiative forces can be used for slowing. As is well 
known from atoms [3, [3] , laser slowing can be effective 
over broad velocity ranges and is insensitive to position, 
so that it can work on a large phase-space volume of 
molecules. Laser slowing can also lead to simultaneous 
longitudinal velocity compression, which is advantageous 
for loading traps. Once trapped, these molecules may be 
further laser-cooled to increase the phase-space density. 
Finally, laser slowing methods appear viable for a variety 
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FIG. 1: Relevant energy levels and transitions in SrF. 

Solid red lines denote slowing lasers, while solid green lines 
denote probe lasers. Dashed lines denote spontaneous de- 
cay channels, and FCFs {fv'v) are labeled if applicable. Blue 
dashed lines indicate decay fluorescence channels used to de- 
termine the beam Doppler profile. 



of species, including some for which current phase-space- 
conserving slowing methods appear ill-suited; this is use- 
ful for addressing the needs of the variety of applications 
envisioned for ultracold molecules (20j . 

Here we experimentally demonstrate deceleration of 
a molecular beam using radiation pressure. This work 
builds upon similar results demonstrating deflection [ll| 
and transverse cooling [l2j of a molecular beam by radia- 
tive force. The crucial enabling feature for radiative slow- 
ing is the ability to scatter > 10* photons without heat- 
ing the internal degrees of freedom of the molecules. Our 
scheme for creating this quasi-cycling transition [lll . [l^ is 
recounted briefly here and depicted in Fig. [T] We employ 
the X2l]+2(u=0, A=1)->A2Hi/2(w'=0, J'=l/2) electronic 
transition of SrF, with t=24 ns lifetime, for cycling and 
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FIG. 2: Schematic of the apparatus. Red (green) lines 
indicate slowing (probe) laser beams. 



slowing. We denote by Wqq 
laser, with wavelength Aqq 



the main cycling and slowing 
and detuning from resonance 



Awqq. The favorable Franck-Condon factors (FCFs) of 
SrF limit vibrational branching j2]| . Separate repump 
lasers, denoted wfg and address residual vibrational 
leakage. Driving an N=\^J'=\/2 transition eliminates 
rotational branching [loj . while a magnetic field remixes 
the resulting dark ground-state Zeeman sublevels [2^ . 
This scheme should allow > 10^ photon scattering cy- 
cles before the bright state population is reduced by 1 /e 



We use the Doppler shift of laser-induced fluorescence 
(LIF) to measure the longitudinal velocity profile (LVP) 
of the molecular beam, as shown in Fig. [21 Our detection 
scheme employs two perpendicular probe lasers, denoted 
Wqq and u^Q, to excite molecules from the X(u=0, 1; 7V=1) 
states (with resolved spin-rotation structure (SRS) and 
hyperfine structure (HFS)) to the A(w=0, J=|) state (un- 
resolved HFS) at a distance L(i = 1350 mm downstream 
from the source. A longitudinal probe laser, denoted 
v\j^, then excites the molecules to the D(z;=0, A^=3, J=|) 
state (unresolved HFS). Monitoring the D— >-X LIF as a 
function of the f^j-, laser frequency yields a Doppler- 
shiftcd LVP free of SRS/HFS, at a wavelength easily fil- 
tered from all laser light. 
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lasers are spatially over- 



lapped countcrpropagating to the molecular beam. To 
address all SRS/HFS levels over a wide velocity range. 
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and lasers have radio-frequency (rf ) side- 



bands with modulation index m = 3.1 and modulation 
frequency 41 MHz < fmod < 44 MHz, unless noted oth- 
erwise. Due to the large frequency extent of the side- 
bands, we do not expect longitudinal velocity compres- 
sion [23[. We note that the dark magnetic sublevels of 
the X(iV = 1) state prevent use of a Zeeman slower. The 
slowing lasers are not chirped [l^l due to the temporal 
extent of the molecular beam pulse (~10 ms). Magnetic 
field coils create an approximately uniform field B=% G 
at an angle Ob =45° relative to the Wqo hnear polarization 
over the length 120 mm ^ L < 1350 mm. A supplemen- 
tary light detector at Lg = 660 mm allows monitoring 
of the LIF from spontaneously emitted photons during 
cycling. 
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FIG. 3: Measured and simulated slowing for different 
detunings of the main slowing laser, a) Measured slowed 
LVP (solid color), control LVP (— ), and velocities correspond- 
ing to the Vqo rf sideband spectrum (gray, with center A). 
The panels are scaled so that all controls have equal heights, 
b) Simulated slowed LVP (solid color), simulated control LVP 
(— ) , and simulated slowed LVP with transverse cooling the en- 
tire length of the slowing (dotted color, scaled by 1/21). The 
gray shaded area indicates the assumed force versus velocity 
profile used in the simulation. The Ai;oo detuning (in MHz) 
is shown in the centered box for each panel set. The scaling 
of the simulated LVP with transverse cooling is chosen to ac- 
curately illustrate how the increased divergence of the slowest 
molecules impacts the apparent loss in our slowed LVP data. 
In particular, the factor 21 is the ratio of the peaks of the 
simulated control LVPs with and without transverse cooling. 
This scaling hence factors out the simple gain in downstream 
flux due to transverse cooling; the different shapes of simu- 
lated slowed LVPs with and without transverse cooling show 
clearly the increased apparent loss of slow molecules in our 
experimental data, merely due to their increased divergence 
(since transverse cooling is not applied in these experiments) 
(Supplementary Discussion 1). 
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Application of the slowing lasers shifts the molecular 
beam LVP, as shown in Fig. [3^ for various Aiiqq. As 
Auqq is tuned towards (w/)/Aqq from the red (where (vf) 
is the mean forward velocity) , the LVP is shifted to lower 
velocities, until Avqq!v{v f) / Xqq] then, when tuned further 
to the blue, the LVP gradually returns to its unperturbed 
state. However, the shift to lower velocities is accompa- 
nied by an increase in apparent molecule loss, which is 
most severe when Avqq^{v f) / Xqq . 

We argue that this apparent loss is due primarily to in- 
creased divergence and transverse heating as the beam is 
slowed. Several other loss mechanisms were ruled out 
as the dominant cause after investigation. Increasing 
the background gas pressure (primarily helium) by 5x 
changed the slowed LVPs little, indicating that back- 
ground gas collisions are not a dominant loss mechanism. 
We investigated possible loss to other rovibrational states 
which could arise from various mechanisms. For exam- 
ple, off-resonant excitation to the A{v = 0,J=^) state 
by the v^q laser or HFS mixing in the A{v ~0, J ~ ^) 
state could transfer population to the dark X(ij=0, A^=3) 
state; stray electric fields could lead to decays from the 
A{v=0, J=i) state to the dark X{v=0; 7V=0, 2) states; or 
the v'^j-, laser could ofii'-resonantly excite molecules from 
the A(?;=0, J=-^) state to the D(w=0,7V=l) state before 
they reach L^- To investigate such mechanisms, we ex- 
plicitly probed the populations of the X(i;=0; N=0, 2, 3) 
and X('!; = l,iV = 0) states and determined that < 10% 
combined total loss could be attributed to such processes. 
Loss to the X(w=3, A^=l) state was not directly measured, 
but was estimated from the observed increase in sponta- 
neous scattering LIF at Ls by adding the V21 repump 
laser; this indicated that molecules cycled through the 
X{v=2, N=l) state ~ 3x before reaching Ld- Together 
with the estimated FCFs [ll|, this yields an estimated 
~6% loss to the X{v=3, N=l) state. Overall, our inabil- 
ity to find evidence of population of dark states, together 
with the observed LVP shift of '^50 m/s, is preliminary 
evidence that our quasi-cycling transition is nearly closed 
for up to ~10^ scattered photons. Additional observa- 
tions provide more quantitative evidence that corrobo- 
rates this conclusion (Supplementary Discussion 2). 

Apparent loss from increased divergence and trans- 
verse heating was modeled via a Monte Carlo simula- 
tion (details in Methods section). Typical simulation 
results, shown in Fig. |31d, indicate that nearly all ap- 
parent molecule loss can be attributed to increased di- 
vergence and transverse heating. According to the sim- 
ulation, the addition of transverse cooling to counter- 
act divergence losses can vastly increase the flux of slow 
molecules. Other observations (e.g., dependence of the 
slowed LVP on laser power for various detunings) cor- 
roborate the basic picture that the observed slowed LVP 
is highly distorted by loss of slow molecules (Supplemen- 
tary Discussion 3). Dependence of the slowing on various 
control parameters was generally weak over small ranges 
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FIG. 4: Slowing with no applied magnetic field. The 

Avqq detuning (in MHz) is shown in the upper right. Oth- 
erwise the representation is the same as in Fig. [S] Note the 
sharp features, in particular the increase in peak height or 
slope of the LVPs; these indicate longitudinal velocity com- 
pression within part of the distribution. This should be con- 
trasted with the smooth LVPs obtained at large B (Fig. [3|. 
Here the Wqo modulation index is m—2.6. 



around their nominal values (Supplementary Discussion 
4). 

Most data was taken with S « 4-9 G and 9b ^ 45°. 
Over this range, the slowed LVPs and LIF at Lg were 
fairly insensitive to the value of B. However, we ob- 
served that Earth's magnetic field. Be, on its own al- 
lows some remixing of the dark Zeeman sublevels. Under 
certain conditions when B = Be (Supplementary Dis- 
cussion 5); qualitatively different behavior was observed; 
namely, sharp features appeared in the LVP, as shown in 
Fig. m Moreover, under these conditions there is clear 
evidence for longitudinal velocity compression: the peak 
and slope of the LVP increases under these conditions for 
certain detunings Awqq. We have been unable to find a 
simple explanation for these features, and full modeling 
of the system (including all ^33 slowing laser frequencies, 
44 molecular sublevels, B-field remixing, coherent dark 
states etc.) is beyond the scope of this paper. How- 
ever, this behavior could potentially be used to compress 
the molecular beam LVP. Ideally this would be done after 
slowing had already removed most of the kinetic energy 
from the beam, e.g., by using an initial region of large 
B for broadband slowing, followed by a second region of 
small B for longitudinal velocity compression and further 
slowing. A slow and nearly monoenergetic beam would 
be ideal for trap loading. 

In summary, we have demonstrated efficient radiation 
pressure slowing of an SrF molecular beam. Under cer- 
tain conditions, we detect ^6% of the initial detected flux 
at velocities <50— . The dominant loss mechanism at 

s 

present is the increased divergence and transverse heat- 
ing of the beam due to the slowing. The addition of 
transverse cooling should provide a much higher flux of 
slow molecules, suitable for loading a trap. It may be 
possible to use a low B-field section to compress the ve- 
locity distribution following the initial slowing. A slow 
molecular beam could be dir ectly loaded into either a 
magneto-optical trap (MOT) [10| or a sufficiently deep 
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conservative trap, using optical pumping as a dissipa- 
tive loading mechanism 25-2^. Furthermore, the pre- 
liminary evidence of little loss during cycling, even after 
>10^ photons have been scattered, invites the possibility 
of moderately long lifetimes for SrF in a MOT. 



METHODS 

We use an ablation-loaded cryogenic buffer gas beam 
source, which provides relatively low initial forward veloc- 
ities, low internal temperatures, and high brightness [2^ . 
[30^ . It produces an SrF molecular beam of 1.2 x 10^^ 
molecules/sr/pulse in the X^E^y2(*'=0i -^=0) state. To mit- 
igate variations of the molecular beam flux, data are taken 
by chopping the slowing lasers on and off between each abla- 
tion shot. The background pressure in the beam propagation 
region is ~2xl0^^ Torr. 

With the foo, fioi smd uli lasers applied, molecules cy- 
cle over the three bright ground states: X(i; = 0, 1, 2; = 1). 
Due to different optical pumping rates into and out of these 
states, the X{v = 0,l; N=l) populations are expected to be 
comparable, while the X(t; = 2, A'' = 1) population should be 
significantly less. We hence employ a scheme to detect pop- 
ulation in both X(u = 0, 1; A'^ = 1) states. These states are 
excited to the A(i;=0, J=3/2) state via the Vqq and w^q probe 
lasers, which are spatially overlapped and intersect the molec- 
ular beam at Ld- Both the Uqq and u^'q lasers have fmod—i^ 
MHz rf sidebands with m = 2.6 to excite all SRS/HFS lev- 
els of the X(u = 0, 1;A^ — 1) states [iJl; since they intersect 
the coUimated molecular beam transversely, they are subject 
to negligible Doppler shift and broadening. The Vgg and u^q 
laser powers are set to drive both transitions with the same 
Rabi frequency. Once excited to the T>{v = 0, N = 3, J — 5/2) 
state by the longitudinally propagating ii^j-, laser, the result- 
ing LIF at Ld, predominantly at 360 nm, is filtered and mea- 
sured by a photon-counting PMT. By scanning the laser 
frequency and noting the Doppler shift (relative to the reso- 
nant frequency of a transverse laser beam), we derive the 
molecular beam LVP. We verified that the detection efficiency 
and measured LVP are independent of whether the molecule 
is detected from X{v^O,N = 1) or X(?; = l,iV = l). Power 
broadening from the laser (~23 MHz FWHM) and mag- 
netic field broadening (~ 18 MHz FWHM) lead to a measured 
broadening of 42 MHz FWHM for the n^^, detection profile, 
resulting in all velocity values having an uncertainty of ±17^ . 

The Vqq, vfo, V21, and v^j~, lasers are spatially overlapped 
using a combination of dichroic mirrors and a polarizing 
beamsplitter (PBS) to produce a single beam with ^ inten- 
sity waist d=3.4 mm (except for the laser with d=4.4 
mm) counterpropagating to the molecular beam. The Vqq, 
Vio, I'll and v^j~i laser powers are 140 mW, 73 mW, 45 mW, 
and 70 mW respectively. The Vqq, vIq, and laser detunings 
from resonance, denoted Awqq, Au^q, and AV21 respectively, 
are first varied iteratively to maximize LIF at Ls- For finer 
tuning, the laser detuning, denoted Aw^^ , is set to reso- 
nantly excite SrF molecules with Vf~50^, and Auoo, Awfo, 
and Aw2i are varied iteratively to maximize the number of 
molecules detected in that Doppler class. Unless explicitly 
noted, Aufo and Ad|i remain at these empirically determined 
values, denoted Av^q^^ and Av^f^ respectively. 

In the Monte Carlo simulation, particles are created at the 



source with randomized velocity distribution matching the 
measured forward and transverse velocity distributions of our 
source [2^. We assume equal detection efficiency over the v^jj 
-ij beam waist at L=Ld- We estimate the force profile using a 
model 5-level system consisting of one excited state and four 
ground states (to match the four SRS/HFS levels). The de- 
generacy of the SRS/HFS levels and the accompanying level 
shifts and remixing within each SRS/HFS level due to the 
applied B-field are not included in the simulation. The satu- 
ration parameter s is calculated for each of the four SRS/HFS 
levels assuming an estimated saturation power of 6 mW/cm^ 
and the known uqo laser rf sideband spectrum. Using r and 
s, classical rate equations are solved to determine the equi- 
librium excited state population fraction, pee, as a function 
of the detuning from the center of the SRS/HFS spectrum, 
A. The dependence of pee is then fit to a Voigt profile. This 
process is repeated for the range of powers dictated by the Vqo 
laser's Gaussian intensity profile. Using the peak values of pee 
for each intensity, we derive an estimate of how the maximum 
scattering rate varies with the distance from the center of the 
slowing beams, r. We finally model the scattering rate R as 
the analytic function 



R(A, 



-Rn 



l+ir/roY 



{wL/2)2+{A-ty 



-dt 



where the normalization A'^i is chosen so that R{Q,Q)=Rmax- 
The parameter values ro = 1.3 mm, a = 3.75, wl = 99 MHz 
and wg = 95 MHz are derived from these fits, without refer- 
ence to the LVP data. The first two parameters control how 
the scattering rate varies with the beam intensity, while the 
latter two characterize the shape of the estimated scattering 
rate as a function of A. Finally, the free parameter Rmax is 
varied manually to fit the LVP data for a variety of A. We 
achieve good agreement with Rmax=2.8 x 10® s~^, consistent 
with our previous observations [ll[|. Additional details of the 
simulation procedure are in Supplementary Methods. 
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SUPPLEMENTARY MATERIAL 
Discussion 1 

The addition of transverse cooling sliifts tlie simulated con- 
trol LVP to lower velocities by < 10^ relative to the unper- 
turbed simulated control LVP. This shift is due to the fact 
that, with transverse cooling, slower molecules are collimated 
closer to the source than are faster molecules, and the likeli- 
hood of their passing through the detection region is there- 
fore greater. Overall, once the difference in scale has been 
accounted for, the shape of the simulated control LVP with 
transverse cooling is similar enough to the shape of the simu- 
lated control LVP without transverse cooling that the former 
is not shown in the main text. The simulated control LVP 
with and without transverse cooling is shown in Supplemen- 
tary Fig. [5] to illustrate this effect. 



Discussion 2 

We use the quantity I^hm, defined as the shift of the half- 
maximum point on the leading edge of the observed slowed 
LVP (versus that of the control LVP), as one simple mea- 
sure to evaluate the effectiveness of our slowing for different 
experimental parameters. Because slowed molecules are less 
likely to be detected (due to increased divergence, etc.), /^hm 
likely provides an underestimate of the actual slowing. With 
Ai)oo=— 260 MHz, providing resonant excitation for molecules 
with w/=i75^, we routinely achieve A£fi\/~45 — 60^. Since 
the SrF recoil velocity is Urc =5.6^^, we interpret this as a 
mean number of photons scattered per molecule {Nsc)~^Q*' , 
roughly an order of magnitude greater than in previous work 

m. 
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Supplementary Figure 5: Impact of transverse cooling 
on simulated control longitudinal velocity profile. The 

simulated control without transverse cooling (black solid line) 
and the simulated control with transverse cooling (black 
dashed line, scaled by 1/21) are very similar; the addition of 
transverse cooling shifts the simulated control LVP by < 10—. 
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Discussion 3 

Other behavior corroborates the basic picture that the 
observed slowed LVP is highly distorted by loss of slow 
molecules. As shown in Supplementary Fig. [5^, LIF at 
Ls peaks at Auqo ~ —160 MHz, whereas /S.hm peaks at 
Auqo ~ —220 MHz. These results are reproduced well by 
the simulation as shown in Supplementary Fig [SJs. We in- 
terpret this discrepancy in peak locations to indicate that the 
greater scattering rate at At>oo = —160 MHz is producing 
slower molecules, but these molecules are more likely than 
faster molecules to diverge before reaching the detection re- 
gion. Further, consider the dependence of Aha/ versus Dqq 
laser power, which is shown for Auoo = —340 MHz, —260 
MHz, and -180 MHz in Supplementary Fig. [7] The data for 
each detuning are fit to a function of the form 



A/ 



Ci 



I + C2P 



where ci and C2 are fit constants and P is the vqq laser 
power. Note that Arm saturates at lower powers when 
Awqq~(i;/)/Aoo, i.e. when the detuning most closely matches 
the Doppler shift of the molecular beam. We interpret this 
to suggest that more slowing is likely occurring for near- 
resonant detuning (~— 180 MHz) than for further red detun- 
ing (~— 260 MHz), but that at higher powers this additional 



Supplementary Figure 6: Dependencies of the half- 
maximum shift and spontaneous scattering fluores- 
cence on laser detuning. Measured (a) and simulated (b) 
dependence of Ahm and LIF at Ls versus Auqo. The LIF 
data is scaled to be the same height as the Ahm peak in 
both panels. 



slowing is not apparent in the data due to increased diver- 
gence and therefore decreased detection probability for the 
slowest molecules. 



Discussion 4 

Overall, the slowed LVP was less sensitive to the ufg and 
V21 laser powers and detunings than to the vqq laser power and 
detuning. At maximum vfo power, Ahm was found to be in- 
sensitive when varying Awfg by ~20 MHz around Awjq°''*. All 
V21 powers >1 mW produced indistinguishable slowed LVPs 
over an even broader range of detunings. Finally, under stan- 
dard conditions where _B = 4— 9 G, the slowed LVP was in- 
sensitive to changes in the ijqo sideband modulation index for 
2.6<m<3.1 and to changes in the dqo sideband modulation 
frequency for 42 MHz < /„od < 44 MHz. 
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Supplementary Figure 7: Dependencies of the half- 
maximum shift on laser power and detuning. Aha/ 
versus dqo laser power for Auoo= -340 MHz (■), -260 MHz 
(•), -180 MHz (a) and associated fits. 



Discussion 5 

In our apparatus, Be ~ .5 G and the angle between Be 
and the Vqq laser polarization is ~ 102°. The data shown 
in Fig. 4 of the main text were taken with a Uqo modulation 
index of m=2.6, where the sharp features were generally more 
pronounced than at the typical m — 3.1 used for most of the 



slowing data. 



Supplementary Methods 

We found it necessary to adjust the size of the simulated 
source to a value larger than the physical cell exit hole diam- 
eter in order to obtain a qualitative match to our data. We 
believe this reflects the fact that collisions outside the source 
increase the effective area from which molecules are emitted. 
Our measurements show the molecular beam transverse ve- 
locity increases significantly (~2x) within ~10 mm from the 
source [2^, an effect we attribute to collisions with residual 
helium in the beam. Similar behavior was noted by a dif- 
ferent group using an apparatus similar to ours [s^. Given 
the measured initial transverse velocity, we calculate the spa- 
tial extent of the molecular beam at the distance where the 
transverse velocity stops increasing (indicating the beginning 
of molecular fiow conditions). We find the calculated spatial 
extent of the molecular beam at this distance to be larger 
than 6.35 mm, the diameter of a hole in a coUimating plate 
34 mm downstream from the source. Ultimately the effective 
source diameter in the simulation was set to 6.35 mm; al- 
though the simulated source remained at the location of the 
physical source exit hole. Using a source diameter of 3 mm 
in the simulation (matching our physical cell exit hole) pro- 
duced apparent losses greater than those measured by a factor 
of two or more for any set of reasonable parameters. 



